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Abstract 

We consider a general formalism to compute inclusive polarised and un- 
polarised cross sections within pQCD and the factorisation scheme, taking 
into account parton intrinsic motion in distribution and fragmentation func- 
tions, as well as in the elementary dynamics. Surprisingly, the intrinsic par- 
tonic motion, with all the correct azimuthal angular dependences, produces 
a strong suppression of the transverse single spin asymmetry arising from 
the Collins mechanism. As a consequence, and in contradiction with ear- 
lier claims, the Collins mechanism is unable to explain the large asymmetries 
found in p ^ tt X at moderate to large Feynman x p . The Sivers effect is 
not suppressed. 



1. Introduction and general formalism 

The inclusive production of large pr particles in the high energy collision of two 
nucleons has been for a long time a crucial testing ground for perturbative QCD; 
in such kinematical regions the partonic degrees of freedom dominate the hadronic 
processes, which can be described in terms of quark and gluon dynamics, coupled 
to non perturbative information - parton distribution (pdf) and fragmentation (ff) 
functions - gathered from other processes and evolved to the proper scale via QCD 
evolution equations. 

In the simplest case this translates into the well known expression: 

^ ,3 = dXadXbdzfa/A{Xa:Q^)fb/B{Xb:Q'^) (l) 

" Pc a,b,c,d-' 



^ ~~~ri — ^fe) <^(s + ^ + It) Dc/c{z, Q^) 
TTz^ at 

= j dXa dXb fa/AiXa, Q^) fb/siXb, Q"^) (2) 

1 da"'^^'^ 

X -J {s,t,U,Xa,Xb) Dc/c{z,Q'^) , 

TTZ at 

which combines all possible elementary QCD interactions ab — > cd, with distribu- 
tion, f{x, Q^), and fragmentation, D{z, Q^), functions: all partonic intrinsic motions 
have been integrated over and the hadrons are considered as composed of collinear 
massless quarks and gluons, each carrying a fraction x of the parent momentum; 
similarly for the final quark fragmentation into a collinear hadron with fraction z 
of the quark momentum. The energy-momentum conservation of the elementary 
interactions, s + i + u = 0, allows to relate Xa,Xb and z, namely, in this collinear 
picture, Xa Xb z s = —Xa t — XbU, where s, i, u (s, t, u) are the Mandelstam variables 
for the partonic (hadronic) process. 

Eq. (1) - taking into account higher order contributions to the elementary in- 
teractions - describes successfully the highest energy cross section data, including 
the most recent ones from RHIC [1]. However, already starting from the pioneer- 
ing work of Feynman, Field and Fox [2], several papers have shown that intrinsic 
transverse momenta fc^'s have to be explicitly introduced into Eq. (1) in order to 
be able to explain data at moderately large p^, for production of pions and photons 
[3, 4]; without them the theoretical (collinear) computations would give results in 
some cases much smaller (up to a factor 10 or even more) than experiment. 

Taking into account intrinsic transverse momenta is not an entirely straight- 
forward matter. In the pure parton model, where partons are regarded as physical 
particles with definite mass (usually assumed to be negligible) , the standard collinear 
parton density fa/A{xa) is simply generalised to fa/A{xa, k±a), where k±a is the par- 
ton momentum perpendicular to the nucleon momentum, and 

fa/A{Xa) = j d'^k±afa/A{Xa,k±a) , (3) 
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where, to be precise, Xa is the hght-cone momentum fraction of parton a inside 
hadron A. Similarly, the fragmentation function is generalised to Dc/c{z, k±c), 
where k±c is the transverse momentum of the observed hadron C with respect to 
the fragmenting parton c. All dynamic partonic calculations are then carried out 
with inclusion of the intrinsic transverse momenta fc^'s. 

This natural generalisation apparently modifies Eq. (1) into: 

d'Pc ^ ' 

^ j dxa dxb dz d^k±a d'^k±b d^k±c S{k±c ■ Pc) fa/A{xa, k±a; Q^) fb/B{xb, k±b; Q^) 



a,b,c,d 

s 



i.2 da"^^"'^ 



-J{k±c) ^ — {s,t,u,Xa,Xb) 5{s + t + u) Dc/c{z,k±c;Q ) , 

TlXaXbZ^S dt ' 

where k^a {k±b) and k±c are respectively the transverse momenta of parton a (b) 
with respect to hadron A (B), and of hadron C with respect to parton c, which in the 
case of light quarks or gluons is taken to be massless. We have formally extended 
our definition of the 2-vector k±c into a 3- vector via the 5-function 5{k±c ■ Pc)- 
Neglecting parton masses, the function J is given by [4] 



[ec + ^/ Pi -klc)' 



Eq. (4) has been widely used in the literature, albeit without the factor J, which 
equals 1 if wc neglect the final hadron mass and k\(j in Eq. (5). Note that the factor 
s/{tiz'^) in Eq. (1) follows from the factor /{ti XaXi,z'^s) in (4) since for collinear 
collisions s = XaXbS. Although it is true that even with fe^, s ~ XaXbS, the use of 
this approximation has been shown by Cahn [5] to lead to azimuthal asymmetries 
which are physically impossible. 

The k± dependent pdf and ff are usually assumed to have simple factorised and 
Gaussian forms, like: 



f,Ux,kr,Q') = Mx,Q')g{k,_) = f^/^{x,Q') (^e-^'''^ , (6) 



so that 



(kl) = , J d'k^ f\/p{x, kr, Q') = Mx, Q') , (7) 

where (3 might depend on x and the energy; it is usually assumed to be flavour 
independent. A similar factorisation is adopted for the k± dependent fragmenta- 
tion functions. The elementary cross sections da/dt depend, via the elementary 
Mandelstam variables s, t and -u, on the intrinsic motions. 

The QCD factorisation theorem implicitly used in Eq. (4) - with unintegrated 
fcx dependent distribution and fragmentation functions - has never been formally 
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proven in general [6], but only for the Drell-Yan process, for the two-particle inclusive 
cross section in e'^e" annihilation [7] and, recently, for SIDIS processes in particular 
kinematical regions [8]. Moreover, in QCD the parton model is a leading-twist 
approximation to the theory, whereas intrinsic transverse effects are of higher-twist 
and should therefore be incorporated in a consistent higher-twist development of 
the theory. Unfortunately, such a treatment is very complicated and introduces a 
whole set of new unknown soft functions and quark-gluon correlations with unclear 
partonic interpretation. 

It turns out, however, that some partonic effects of transverse momentum are 
surprisingly large and can generate phenomena which would be impossible to repro- 
duce in the coUinear treatment: 

• the presence of an intrinsic k± alters the relationship between the light-cone 
momentum fraction x of the parton and the Bjorken x^j^ so that X ^ XBj- 
Although the shift is small and proportional to k']_/(x^/s)'^, it can have a 
substantial effect in the region of x where the parton densities are varying 
rapidly. This is a kind of enhanced higher-twist effect and can lead up to an 
order of magnitude change in a cross section. Similarly, due to intrinsic motion, 
the partonic scattering angle in the pp cm. frame might be much smaller 
than the hadronic production angle, thus enhancing the large pt inclusive 
production of particles. 

• In the presence of transverse momentum, certain spin-dependent effects can be 
generated by soft mechanisms and can be used to understand the large trans- 
verse single spin asymmetries (SSA) found in many reactions like -\- B ^ 
C-\-X and the large hyperon polarisations in processes hke A-\- B ^ -\- X . 
At leading twist there are 4 such soft mechanisms, often referred to as "odd 
under naive time reversal": 

a) Sivers distribution function [9]: in a transversely polarised nucleon with 
momentum p and polarisation vector P, the number density of quarks with 
momentum (xp, k±) is allowed to depend upon P • {p x k±) ; in other words, 
the Sivers distribution function represents the azimuthal dependence (around 
p) of the number density of unpolarised quarks inside a transversely polarised 
proton. 

b) Collins fragmentation function [6] : in the fragmentation of a transversely 
polarised quark with momentum p^ and polarisation vector Pq, q C-\-X, the 
number density of hadrons C with momentum {zp^, k±) is allowed to depend 
on Pq ■ {pg X k±); in other words the Collins fragmentation function represents 
the azimuthal dependence (around Pq) of the number density of unpolarised 
hadrons resulting from the fragmentation of a transversely polarised quark. 

c) Boer-Mulders distribution function [10] : in an unpolarised nucleon a quark 
with momentum [xp, k±) is allowed to have a non-zero polarisation along 
p X k±; that is, the Boer-Mulders distribution function represents the az- 
imuthal dependence (around p) of the number density of transversely polarised 
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quarks inside an unpolarised proton. 

d) polarising fragmentation function [11, 12]: in the fragmentation of an unpo- 
larised quark with momentum a final spin 1/2 hadron C with momentum 
{zpq,k±) is allowed to have a non-zero polarisation along x k±; that is, 
the polarising fragmentation function represents the azimuthal dependence 
(around p^) of the number density of transversely polarised hadrons resulting 
from the fragmentation of an unpolarised quark. 

It should be noted that in the pure parton model, where partons are treated as 
physical free particles, all these effects vanish [13] . 

In the present paper we study transverse single spin asymmetries, in p ^ n X 
processes, taking into account all parton intrinsic motions in initial and final hadrons 
and in the elementary dynamics. This generalises previous work in which only the 
one fe_L essential to the mechanism was taken into account, either in the initial 
polarised nucleon (Sivers effect) or in the final quark fragmentation (Collins effect), 
and the k± distribution was somewhat simplified into essentially a two-dimensional 
5- function [14, 15, 16]. 

For the reasons explained above, we have not attempted to construct a fully 
consistent next-to- leading- twist treatment. Our strategy is to keep only the en- 
hanced higher-twist terms and to calculate partonic helicity amplitudes as if the 
partons were particles. We believe this approach is physically meaningful since it 
takes into account the most important higher-twist terms in the cross section and 
the asymmetry. Three of the above spin effects, a)-c), can contribute to pion SSA, 
but in this paper we wish to explore the generation of SSA due to the existence of 
the Collins fragmentation function alone, for which there is some evidence in the 
polarised lepto-production data of the HERMES collaboration [17, 18]. The Bocr- 
Muldcrs effect can also contribute to transverse single spin asymmetries but, at least 
forp^p — > ttX processes, it would contribute mainly at negative xp values, whereas 
data are in the positive xp region. The Sivers effect is also relevant, and has been 
studied in a parallel paper [4]. In fact we shall show that the consistent treatment 
of all intrinsic partonic motions induces a major suppression of the contribution to 
the asymmetry due to the Collins mechanism and renders it incapable of producing, 
by itself, the kind of asymmetries measured in p^ p ^ n X reactions [19, 20]. 

This result modifies the conclusions of Ref. [15, 16], where the Collins contribu- 
tion to SSA in p — > TT X processes was computed adopting a simplified kinematical 
configuration: it appeared that the Collins fragmentation function could, although 
with some difficulty, explain the E704 data [19]. Note that the results on SSA ob- 
tained using the Sivers distribution function, with a similar simplified kinematical 
configuration [14] are, instead, essentially confirmed by the exact treatment of all 
intrinsic partonic motions [4]. 

In order to study spin asymmetries we have to introduce spins in the QCD hard 
scattering processes. Eq. (1) holds also for polarised processes, {A, Sa) + {B, Sb) 
C + X [21], provided one introduces in the factorisation scheme, in addition to the 
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distribution functions, the helicity density matrices which describe the parton spin 
states. This can be done also for Eq. (4) with the result: 



E^^^{A,Sa)HB,Sb)^C+X ^ dXgdXbdz 



J {k±c) pI[X^ fa/A,SA {Xa, k±a) p'lx" -^VB-Sb i^b: k±b) (8) 

Mk,x„K,K <,a,;A,,a^ S{s + i + u) bl±{z, k^c) , 

where we have used the notation {A} to imply a sum over all helicity indices. In 
Eq. (8) Px^w^ is the helicity density matrix of parton a inside the polarised hadron 
A whose polarisation state is generically labelled by Sa (for spin 1/2 particles this 
means longitudinal or transverse polarisation); similarly for parton h inside hadron 
B with spin Sb- The Ma^,a^;A„,a^'s are the helicity amplitudes for the elementary 
process ah cd, normalised so that the unpolarised cross section, for a coUinear 
collision, is given by 

- 16^4, , ■ 

'~ A A' 

■^\^',x''^ {^1 k±c) is the product oi fragmentation amplitudes ior the c — >• C+X process 

D^o^^^ =4 Vx X -xVl w .y , (10) 

where the X stands for a spin sum and phase space integration over all unde- 

JX,Xx 

tected particles, considered as a system X. The usual unpolarised fragmentation 
function Dc/dz), i.e. the number density of hadrons C resulting from the fragmen- 
tation of an unpolarised parton c and carrying a light-cone momentum fraction z, 
is given by 



Dc/c^ 



= ^ E jd'k^cDl±{z.k^c). (11) 



Eq. (8) can be formally simplified, showing its physical meaning, by noticing 
that: 

A,Sa B,Sb 



PkA Pkk ^a..a,;A„,a, Ml,^,.,,^,, ^ (c) = PA..A, Trp'(c) , (12) 



where Pa^,a^ is the normalised helicity density matrix of parton c produced in the 
ab — s> cd process, with initially polarised partons a and 6; the normalisation factor 
Trp'(c) is related to the polarised cross section for a coUinear coUision: 

^2 - (a,So)+(b,Si,)^c+a! 

IVp'(c) = (327r^5^) , (13) 

dt dcp 
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where is the azimuthal angle of parton c in the partonic center of mass frame. 
Moreover, 

E PK,K d'Ck^'^ ^^o) = Dc/csA^, k^c) , (14) 

is just the fragmentation function of a polarised parton c, with spin configuration 
Sc, into a hadron C, whose spin is not observed. 
Using Eqs. (12)- (14), Eq. (8) can be written as: 



d^Pc 



J2 j dxa dxb dz d'^k±a d'^k±b d^k_ic S{k±c ■ Pc) 



C a,b,c,d 

fa/A,SAi^a, k±a) fb/B^Ssi^b, k±b) ^ J{k±c) 

Xa-^b^ ^ 

^2^{a,Sa)+{b,si,)^c+d 

-pr-. 5{s + t + u)Dc/c,sS^,k^c) , (15) 

dt d(p 

which is the analogue of Eq. (4) in the polarised case. 

Eq. (15) shows clearly the factorised structure and the partonic interpretation: 
inside polarised hadrons one has polarised partons with spin configurations and 
Shi which interact via pQCD processes, leading to a final polarised parton, with spin 
configuration Sc, which fragments into the observed final hadron. For the initial and 
final step - the determination of the parton polarisation from the hadron polarisation 
and the fragmentation of the polarised parton - one has to rely on distribution and 
fragmentation functions; some of them are known from other processes or from the- 
oretical models and some of them, in particular when allowing for intrinsic motions, 
are new and unexplored. 

Although Eq. (15) has a simple physical interpretation, it is more convenient to 
study the scattering process with the helicity formalism of Eq. (8); when dealing 
with helicities and helicity density matrices all spins have a well defined interpreta- 
tion concerning their directions [22], and this is crucial if we arc taking into account 
all parton transverse motions, so that there are several transverse spin directions. 
Since the direction of motion of the parton does not coincide with that of its parent 
hadron, the longitudinal and transverse direction of the parton spin will also be 
different from the longitudinal and transverse direction of the parent hadron spin. 

The partonic distribution is usually regarded, at Leading Order, as the inclusive 
cross section for the process A — > a + X; therefore the helicity density matrix of 
parton a inside a hadron A with polarisation Sa can be written as 
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having defined 



^A^.a:, = f ^^_^^>A„,A^^;A^ -^A^.A^^iA:, > U«J 

and where the x stands for a spin sum and phase space integration over all 



undetected remnants of hadron A, considered as a system Xa and the ^'s are the 
helicity distribution amplitudes for the A ^ a + X process. 

Notice that Eq. (17) relates the helicity density matrix of parton a to the helicity 
density matrix of hadron A. The helicity density matrix describes the spin orienta- 
tion of a particle in its helicity rest frame [22]; for a spin 1/2 particle, Tr {(Tip) = Pi is 
the i-component of the polarisation vector P in the helicity rest frame of the parti- 
cle. In this sense Eq. (17) relates the hadron polarisation to the parton polarisation, 
which have both to be defined and interpreted in the proper rest frames. 

The distribution function of parton a inside the polarised hadron A,Sa is given 

by 



fa/A,SA{^a,k^a)= E PaS'/a'Ia'. (19) 



Aa.A^,A'^ 



and the usual unpolarised distribution function fa/A{xa), i.e. the number density of 
partons a inside an unpolarised parton A, carrying a light-cone momentum fraction 
Xa, is given by 

fa/A{Xa) = l E jd'k^<^^i:i\- (20) 
AoiA^ 

Similar expressions for the fragmentation process have already been introduced 
in Eqs. (10) and (11). 

By using Eq. (17), Eq. (8) can be written as 

Ec da'^^'^^'^+''^'^''^^'^+^ r dxadx^dz 

Ak^c) ptX p"'X ^a^aI <,a,;A,,a^ 5{s + i + u) dI±. (21) 

Eq. (21) contains all possible combinations of different distribution and fragmenta- 
tion amplitudes: these combinations have partonic interpretations and are related 
to the k± and spin dependent fragmentation and distribution functions discussed 
above and, for example, in Refs. [23] and [24]. Notice that, even though Eq. (4), for 
the unpolarised cross section, looks intuitively correct and convincing, Eq. (21), if 
Collins and Boer-Mulders effects arc operative, will yield a different result, i.e. with 
p{ y = {l/2)6x^^x'j {I = A,B), Eq. (21) contains terms not included in Eq. (4), 
that is the terms off-diagonal in the parton helicities. We have checked numerically 
that these contributions are negligible in the unpolarised cross section. All this will 
be discussed in detail in a forthcoming paper [25] , where all contributions to single 
and double spin asymmetries will be examined, together with the parity and k± 
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properties of the distribution and fragmentation amplitudes. Here we are only con- 
sidering the process p — > tt X and are focussing only on the contribution of Collins 
mechanism [6] , that is the azimuthal dependence of the number of pions created in 
the fragmentation of a transversely polarised quark. The unpolarised cross section 
will be computed according to Eq. (4), taking into account the intrinsic transverse 
motion of all partons (see also Ref. [4]). 

2. Single Spin Asymmetries and Collins mechanism for pion production 



Let us then consider the processes p\p^)p — > ttX; we study them in the pp 
center of mass frame, with the polarised beam moving along the positive Z-axis and 
the pion produced in the XZ plane with {p„)x > values. The t (J,) is defined 
as the +Y{—Y) direction. We then have, with Sa = T) i) with an unpolarised 
hadron B {Sb = 0), 

p^'^' - 1 M P^'' --(^ M • (22) 

P^a^'a ~ 2\±i 1 ) ^^b^'b ~ 2 1^ 1 ) ^"^"^^ 

The computation of the single spin asymmetry 

da^ — da^ _ 

requires evaluation and integration, for each elementary process ah — > cd, of the 
quantity [see Eq. (21)] 

E(5a, Sb) = pifX ^>^aX P^b% ^^bX ^K^a-^a,\ Dl^,K ' (24) 

{A} 

where £>J a' defined as in Eq. (10), for pion production. Prom Eqs. (22) and 
(24) one has that the numerator of Ajv is proportional to 



E(T,0)-E(i,0) = X 



Ma A -A A Ml . . w w Dl y , 



(25) 

while the denominator contains: 

E(T,0) + S(i,0) = E \ [Fl:i''^ + F'_f^] ^Atl^W,A,<,A,;A,,A^4:,A,- (26) 

In the equations above and in the sequel + and — stand for +1/2 and —1/2 helicities, 
when referring to nucleons or quarks, and for +1 and —1 helicities, when referring 
to gluons. 

As we have said, in this paper we are focussing solely on the Collins mechanism 
and we do not consider all possible contributions to /l^v, which will be discussed 
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elsewhere [25]. Therefore, we do not consider the possibihty of finding transversely 
polarised quarks inside the unpolarised proton B [10] or the possibility of having 
different total numbers of quarks, at different k±_ values, inside the transversely 
polarised proton A [9]. This does not imply that these other effects which are 
negligible for the impolariscd cross section are negligible for the SSA; simply that 
we wish to explore to what extent the Collins mechanism alone is able to explain the 
measured transverse single spin asymmetries. As a consequence, the F-terms off- 
diagonal in Af,, (while diagonal in A^, A^) and the F-terms off-diagonal in A^, A^ 
(while diagonal in A^, A^) will be neglected. The Collins mechanism corresponds to 
the terms off-diagonal in the fragmenting quark helicities A^, A^. Taking all this into 
account, a partial summation in Eq. (25) obtains 



s(T,o)-s(i,o) = E^^{ 



FZ'Z - f: 



a/A 

where we have exploited the fact that, by parity invariance, 

+ = (28) 

independently of the value of A^. 

The same procedure applied to Eq. (26) reveals that the denominator ol is 
just twice the unpolarised cross section, as given by Eq. (4). 

Eq. (27) can be further simplified by exploiting the dynamical and the parity 
properties of the helicity amplitudes appearing in it. This requires some careful 
considerations. 



• Whereas the hadronic process p ^ t: X takes place, according to our choice, 
in the XZ plane, all other elementary processes involved: A{B) a{h) + X, 
ah —>■ cd and c ^ it + X, do not; all parton and hadron momenta, p^, Pq, 
have transverse components k^a, k±b, k^c and this complicates remarkably 
the kinematics. For example, the elementary QCD process ah — > cd, whose 
helicity amplitudes are well known in the ah center of mass frame, is not, 
in general, a planar process anymore when observed from the pp center of 
mass frame. Similarly, as we commented, the spin properties described by 
helicity density matrices have clear physical interpretations in each particle's 
own helicity rest frame, but not necessarily in the pp center of mass frame. 
Of course, one can always boost and rotate from one frame into another, but 
this introduces phases in the helicity amplitudes, which have to be properly 
accounted for. 
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We refer all angles to the pp cm. frame, in which = {9i, 4>i), {i = a, b, c, d). 
Then the distribution functions of the polarised proton A describe processes 
taking place in the plane defined by Z and the p„ direction, {9a, (t>a)- Therefore 
[22, 25] 

^K,\^^;X^{Xa, k±a) = ^X,,X^^;X^{Xa, k±a) exp[iA^0j (29) 

and 

F^^-y, i^a, k^a) = {xa, k^a) exp[z(A^ - X'M , (30) 

A A' ^ A A' 

where k_i_a — \k±a\', "\? (xa, k±a) has the same definition as F. (xa, k±a), 
Eq. (18), with replaced by JF. 

The parity properties of ^\^,x^ ;A^ (a^a, ^±a) are the usual ones valid for helicity 
amplitudes in the 0a = plane [22], 



^-A„,-A,^;-A, = V (-1)^— (_1)A.-A.+A,^ ^A„,A,^;A, , (31) 

where 7] is an intrinsic parity factor such that rf = 1. These imply: 

p-K-K ^ (_i)2(5a-s„) (_i\{Xj^-x^)+{x'^-x'^) pK,K C32) 

Let us consider now the elementary partonic amplitudes. As already remarked, 
the hard partonic interactions, a[p^ + 6(^5) <^'{Pc) + d{Pd)j take place out 
of the XZ plane, which we have chosen as the plane of the overall p ^ n X 
process. One could compute the helicity amplitudes for these generic processes 
among massless particles using techniques well known in the literature, like 
those explained in Chapter 10 of Rcf. [22]. On the other hand, the explicit 
expressions and the parity properties of the helicity amplitudes M^, which 
apply when the elementary scatterings occur in the ab cm. frame, in the XZ 
plane, are well known. Therefore, rather than computing directly the generic 
helicity amplitudes M, we prefer to relate them to the known amplitudes M°. 

To reach the simple configuration of the M° amplitudes, starting from the 
generic configuration p^, Pi^, we have to perform a boost in the direction de- 
termined by (p^ + pj) so that the boosted three-vector (p^ -|- p[) is equal to 
zero. This will provide us with a c.m.-like reference frame S' where the partons 
a and b collide head-on. Here the parton a and the parton c, resulting from 
the hard interaction between a and b, will have directions identified by {9'^, 0^) 
and {9'^, 0^) respectively. In general, the parton momenta in S' are related to 
the initial ones (before the boost) by: 




P^ = P^-^T-T^[-V^+P^ (^3) 



where i = a,b,c,d and = (g°, q) = Pa + Pb- 
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We need now to perform two subsequent rotations, one around the Z axis by 
an angle and one around the Y axis, by an angle 6'^, such that the collision 
axis of the two colliding initial partons turns out to be aligned with the Z axis. 
We call this frame S". 

Under these boost and rotations the helicity states and consequently the scat- 
tering amplitudes acquire phases, Ca,b,c,d and ia,b,c,d' 

My X , , = M^'\ Q-i{>'a^a+\£.b-Kic-\a.^d) Q-i[{>'a-\)la-{K->'d)^'=\ ('34) 

where and (j = a, 6, c, d) are defined by [22] 

cos 6q sin 6j — sin 6q cos 9j cos{(f)q — (pj) . . 

^ ^^^^ 

sine, = ^^^h^i^l^, (36) 
sin Oqp. 

ij = V- + e ■ , (37) 

cos 9q sin 9', — sin 9q cos 9', cos(0g — 0' ) 
= .ina,,: ^ 

sine = -^'""^'y.-^i). (39) 

■' sm aigp^. 



and 
where 



COS^^-COS^^COS^pj^p. 

cos 77,- = —TIT— 7^ (40) 



sin 6*^ sin 9„' , 



. , sing;^ sin(0; -(/>;.) 

smTy . = — , 41) 

sm 9w J,' 

and the polar angles {9j,(j)'j) arc determined via Eq. (33). Here 9p^p. (0 < 
dpiPj ^ tt) is the angle between p,j and p^-, and so on. Notice that rj'^ = 0. 

In the S" frame the direction of the parton c is characterised by an azimuthal 
angle (f)" given by 

, ,// ^ sing^ sm{(p[ - (f)'^) 

'^^ sin^;^cos(0;^-0;,)cos^;,-cos^;^sin^;, ' ^ ' 

A final rotation around Z of an angle 0^' will then finally bring us to the 
canonical configuration in which the partonic process is a cm. one in the XZ 
plane. This introduces another phase. As a result of the performed boost and 
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rotations the elementary scattering amplitudes computed in the hadronic cm. 
system (the one where we are studying the hadronic cross section) are related 
to the helicity amplitudes computed in the partonic cm. system (in the XZ 
plane, 0c = 0) by: 



g-i(A,^a+A,€i,-A,Cc-A^Cd) g-i[(A„-Aj)ea-(A,-AJ^c] ^i{K-\)4-'L 



(43) 

with 0^', and defined in Eqs. (35)-(42); Eq. (33) allows to fully express 
the amplitudes in terms of the -p'p cm. variables Pj. The parity properties of 
the canonical cm. amplitudes M° are the usual ones: 



A,,-A,;-A„,-A, = ^aWc% 



where r]i is the intrinsic parity factor for particle i. 



(44) 



Let us finally consider the fragmentation process. We take as independent 
variables, in the pp cm. frame, the four-momentum of the final hadron 



Vc 



= {Jpt + Pl,Pt,0,Pl) (whose three-momentum, according to our 



choice, lies in the hadronic XZ plane and where we neglect the pion mass), 
the intrinsic transverse momentum k±c = ^^tt = {k±n,0k^^,4>k^^) of the final 
pion with respect to (fcj_7r ■p^ = 0), and the light-cone momentum fraction 

The parity properties of the fragmentation amplitudes, Eq. (10), are simple - 
analogous to the ones for the distribution amphtudes, Eqs. (31) and (32) - in 
a frame 5'''^ in which the parton c moves along the Z'^-axis. This frame can be 
reached from the hadronic pp frame by performing two rotations: first around 
Z by an angle (pc and then around the new F-axis by an angle 6c, which brings 
the 3-momentum of parton c along the new Z^-axis. In the frame the 
azimuthal angle (f)^ identifying the direction of the final detected pion (which 
coincides with the azimuthal angle of k±T^ in S^) is given, in terms of our 
chosen pp cm. variables, by 



tan(/)f = ± 



Pt 




(45) 



where E^^ = \Jpt + Pl is the energy of the final pion. 

The analogue of Eqs. (29), (30) and (32), for the fragmentation of a parton c 
into a pion, reads 



'Dx-xS^^kx^) exp[iAc0f] , 



^LAi = ^LAiexp[z(\-A'J0^], 



(46) 
(47) 
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with the parity relationships 



(48) 



where DJ^ y is defined according to Eq. (10), in the case in which the hadron 
C is a spinless particle (pion), 



(49) 



By exploiting the above angular and parity relations, Eqs. (30), (32), (43), (44), 
(47) and (48), we can now further simplify Eq. (27). One obtains: 



E(T,0)-E(i,0) = -^E[A/B< 



0* 



{>-} 



{F+r COs[<f>a + € - - la + 2 X^^, + I + )] (50) 
- F+- COS[0„ + + la- 2 A,(ee + lc + 0f)]} , 

where we have also used the fact that partons a and c, carrying transverse polarisa- 
tion, are quarks or antiquarks, that is = •^c = 1/2- 

Let us finally perform the remaining sum over helicities in Eq. (50). The only 
types of elementary interactions contributing are g^gt — > qcQd (generically denoted 
as qq) and qg — > qg (generically denoted as qg)^ where = u, d, s, u, d, s and so on. 
The only independent helicity amplitudes M° for the qq processes are: 



(51) 



and, for the qg processes. 



i)q9 





+,-1; 



(52) 



At Leading Order all such amplitudes are real. 

On summing over {A} Eq. (50) gives, for qq processes. 



[E(T,0)-E(i,0)^ = {F+:{Xa,k^a)cOs[(Pa + <P:-Ca-L + Cc + L + <P^] 

- F+-(Xa, k^a) COS[0„ - ct>l + ^a + la - - L - 0f ]} (53) 

z;k±a, k±b, fe±7r) -2iDl_ {z, k±^) 

and, for qg processes 

[E(T,0)-E(i,0)]^^ = {F+l{Xa,k^a)cOs[<f>a + <f>'i-^a-L + ^c + L + €] 

- F+;(Xa, k^a) COS[0„ - (/.'J + + la - - Ic - '/'f ] } (54) 



fe_La, fc_L6, fc_L^) -2iD!^_(2;, /C. 
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The product of amplitudes appearing in Eqs. (53) and (54) are given by: 



M^M^ = ^^^^ (q^q^^^q^qs) (55) 

where a, /3, 7 and 5 are flavour indices. Notice that in the above expressions aU 
the dependences on the angles in the distribution and fragmentation functions are 
explicit and the functions F , f and D do not depend on angles any more; the 
elementary amplitudes depend on angles via the Mandclstam variables s, t and u. 
Notice also that the qq and qg contributions have exactly the same structure, the 
difference l)eiug only in the parton b distribution and in the elementary processes. 

From Eqs. (21), (24), (53) and (54) the numerator of the single spin asymmetry 
Atv, under the assumption that only Collins effect contributes, is given by {b,d can 
be either quarks or gluons): 



(56) 



daP^ P^'' ^ daP^ p^"" ^ _ 

d^P^ d^p^ 

I dx dx.dz ^3^^^ 

X J{k_i_^) 6{s + i + u) 

X {F+r (X„, k^a) COs[(/)„ + -ia-L + ic + L + (t>"] 
- F+-{Xa, k^a) COs[(/)„ - + + L " " L " 0f ] } 

fc_La, fc±6, fe±7r) , , \-2iDl_{z,kx_^)\ . 

L J qaO^qcd L J 

A few comments are in order. 

• All angles appearing in Eq. (56) can be expressed in terms of the pp cm. 
integration variables, via Eqs. (33), (35)-(42) and (45). 

• Prom Eqs. (48) and (49) one can see that is a purely imaginary quantity. 
The Collins fragmentation function [6, 15, 24, 27] 

-2iDl_ = 2 lmDl_ = A^D^/^, , (57) 

has a simple interpretation in the frame in which the quark moves along the 
Z direction, with spin parallel (q^) or antiparallel (q^) to the F-axis, while 
the q ^ TT X process occurs in the XZ plane: it gives the difference between 
the number density of pions resulting from the fragmentation of a quark q^ 
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and a quark qK In the pp cm. frame the quark transverse spin direction is 
not, in general, orthogonal to the n X plane and this reflects into the 0^ 
dependence appearing in Eq. (56). 

The product of elementary amplitudes Mf'Mg, see Eqs. (51) and (52), is, in 
a frame in which the partonic cm. scattering plane is XZ, simply related to 
the spin transfer cross section: 



1 



IGtts^ 



M°M° = . . 58 

^ dt dt ^ ^ 



Again, the parton intrinsic motions give, in general, more complicated, non 
planar conflgurations for the elementary scatterings, which induce dependences 
on the angles ^j, ^j, rjj and (f)'^. 

• The distribution terms F^z{xa, k±a) and F^^{xa, k±a) are related to the dis- 
tribution of transversely polarised quarks inside a transversely polarised pro- 
ton; these transverse directions can be different for protons and quarks [25]. 
Without any intrinsic motion, only the F^z{xa) distribution would be present, 
coinciding with the transversity distribution hi{xa) [24]. 

• Note that if one takes into account intrinsic motions only in the fragmentation 
process, assumed to occur in the XZ plane [fc_La = k_n, = 0, (fc_L7r)y = 0, which 
implies all phases to be zero], one recovers the expression for the numerator 
of An (aside from the factor J) used in Refs. [15, 16]. 

We can now use Eqs. (56) and (4) to compute the SSA An = {da^ — da^)/2 da. 



3. Attempts to fit the data: suppression of the Collins mechanism 

As noted earlier, it was previously believed that the remarkably large SSA found 
e.g. in the E704 experiment [19] could be generated by either the Sivers [14] or 
the Collins mechanisms [15, 16]. However, to avoid handling the very complex 
kinematics and having to deal numerically with 8-dimensional integrals, only the 
one essential intrinsic responsible for the asymmetry, was taken into account 
in these studies. We now believe that the phases involved, when the kinematics is 
treated carefully, are crucial, and, as we shall see, lead to a large suppression of the 
asymmetry due to the Collins mechanism. As explained in [4] there is little or no 
suppression of the asymmetry due to the Sivers mechanism. 

In order to demonstrate the extent of the suppression we shall choose for the 
unmeasured soft functions in Eq. (56) their known upper bounds. Let us first write 
these functions with the notations of Refs. [11] and [24] (details will be given in 
[25]): 

F^l(x, k±) = hi(x, k±) = hirix, k±) + — ^ h^ri^, k±) (59) 
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F^+{x,k^) = ^hiAx,k^) (60) 

2k 

-2iDl_{zM) = A''D,,^,{zM) = ^^Hi\zM). (61) 

where Mp and M^- are respectively the proton and pion mass. The following posi- 
tivity bounds hold [28, 29]: 

\hi{x, /c_l) I < ^ [q{x, k±) + Aq{x, k±)] = q+{x, k±) (62) 
:^^\hiT{x,k±)\ < -[q{x,k±) - Aq{x,k±)] = q_{x,k±) (63) 

|A^L>,/,t(^,A;^)| < 2D,,,{zM)- (64) 

In our numerical estimates we adopt for all the unmeasured soft functions the 
above maximum possible values, and, moreover, adjust their signs so that the con- 
tributions from the valence flavours (up and down) reinforce each other in the tt"*" 
reaction, producing a maximally large positive A'^j^ . By isospin invariance it then 
turns out that this choice also produces a maximally large negative A^^ . To be 
precise, we have computed the SSA, Ajsi = {da^ — da^)/2da, via Eqs. (56) and (4), 
with the following choices: 

• For the transversity pdf F^z{x, k±) = hi{x, k±) and its companion h^j, we have 
only considered up and down quark flavours, without any sea contribution. We 
have saturated Eqs. (62) and (63): 

hi{x, k±) — u+{x, k±) hf{x, k±) — —d+{x, k±) (65) 

^ ' hi^ix, k^) = -u4x, k^) hiSix, k^) = +d4x, ks.) . (66) 



2M2 ^ ' -V , 



One naturally expects, for valence quarks, positive values for /i" ^-nd negative 
ones for hf; the relative signs between hi and hj-j, are chosen in order to max- 
imise the sum of their contributions in Eq. (56). The x and kj_ dependences in 
the unpolarised and polarised pdf arc factorised assuming the same Gaussian 



form as in Eq. (6), with y (A;^) = 0.8 GcV/c [4]. For the x-dependence of the 
unpolarised pdf we have adopted the MRSTOl set [30] and for the polarised 
pdf either the LSSOl set [31] or the LSS-BBS set [32], as two examples of very 
different choices. We have used the same QCD evolution scale as in Ref. [4] . 

We have chosen the z-dependence of the Collins function in such a way as 
to maximise the effects. Let us consider the production of 7r"'"'s: since the 
dominant partonic contribution at large Xp is ug ^ ug, for which the product 
of elementary amplitudes Mj^M^ is negative, see Eqs. (55), in order to get a 
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positive An we need a negative w-quark Collins function. That is, we satisfy 
the positivity bound (64) with: 

A^D,+/„T {z, k^) = -2D,+,^{z, k^) . (67) 

We consider here also the contribution of the sub-leading channel dg —>■ dg 
(neglected in Refs. [15, 16]); as it enters with a negative hf, in order to add 
all contributions, we use for the non-leading Collins function 

A^D,+/,T {z, k^) = +2D^y,{z, k^) . (68) 

In this way also for 7r~'s is maximised in size (by isospin invariance). 
For 7r°'s we take, exploiting isospin symmetry, 

A^L>,o/,T = \ (A^i?.+/«T + ^''D^+/d^) = \ {-2D,y^ + 2D,+,a) , (69) 

where q = u,u,d,d and which still fulfills the bound (64). The z and k± de- 
pendences of the unpolarised fragmentation functions are also factorised, with 
the same Gaussian dependence as in Rcf. [4], which introduces a 2;-dependent 
(A;^) value, smaller than the constant (A:^) value assumed for the pdf. This 
value allows a good understanding of the unpolarised cross sections; we have 
explicitly checked that increasing the ff {k^) does not change significantly our 
present results (while spoiling the agreement with the unpolarised cross sec- 
tions). The ^-dependent unpolarised ff are taken cither from Kretzer [33] or 
from KKP [34, 4] , as typical examples of two different sets. 

With the above choices, Eqs. (59)-(69), we can (over) estimate the maximum 
value that, within our approach, the Collins mechanism alone contributes to the 
SSA in p ^ 71 X processes. The results are presented in the four plots of Fig. 1, 
which show (^Ar)^^™'^ as a function of xp, at pt = 1-5 GeV/c and \/s ~ 19.4 
GeV: this is the E704 kinematical region and a comparison with their data [19] 
is shown. The only difference between the plots is given by different choices of 
the polarised distribution functions and/or the unpolarised fragmentation functions. 
Four different combinations are possible: two different sets of polarised pdf, LSSOl 
[31] or LSS-BBS [32], and two different sets of unpolarised ff, Kretzer [33], or KKP 
[34]. The four combinations exhaust all possible features of choices available in 
the hterature. The results clearly show that the CoUins mechanism alone, even 
maximising all its effects, cannot explain the observed SSA values; its contribution, 
when all proper phases are taken into account, fails to explain the large E704 values 
observed for AJ^ and A'lf at large xp- 
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4. Comments and conclusions 

Wc have developed a consistent formalism to describe, within pQCD and a fac- 
torisation scheme, the inclusive production of particles in hadronic high energy col- 
lisions; all intrinsic motions of partons in hadrons and of hadrons in fragmenting 
partons, are properly taken into account. Such a scheme has been applied, in a 
parallel paper [4], to the description of several unpolarised cross sections and to 
the computation of SSA in p ^ ir X processes, generated by the Sivcrs mecha- 
nism alone. In this paper we have again considered SSA in p^ p — > tt X processes, 
but focussing on the contribution of the Collins mechanism alone. Previous work 
[14, 15, 16], performed in a similar scheme with simplified kinematics, showed that 
both the Collins and the Sivers mechanisms, could alone explain the observed data 
on SSA. 

Such a conclusion has now to be modified: while properly chosen Sivers distri- 
bution functions could still explain the data [4] , there are no Collins fragmentation 
functions able to do that, as Fig. 1 shows. The failure of the Collins mechanism, 
when all partonic motions are included, can be understood from the complicated 
azimuthal angle dependencies in Eq. (56): the many phases arising in polarised 
distribution and fragmentation functions, and in polarised non planar elementary 
dynamics conspire, when integrated, to strongly suppress the final result. 

The situation with the Sivers contribution alone is much simpler, as the partons 
participating in the elementary dynamics and in the fragmentation process are not 
polarised. As a consequence, the phase structure of the numerator of An, in Sivers 
case, contains only one phase, the Sivers angle (see Eqs. (44) and (45) of Ref. [4]). 
Its integration, coupled with the dependence of the elementary dynamics on the 
same angle, does not significantly suppress the result. In this case, the simplified 
kinematics of Ref. [14] contains the main physical features of the mechanism and 
gives a reasonably accurate computation of A^. 

Our results show, once more, the importance and subtleties of spin effects; all 
phases have to be properly considered and they often play crucial and unexpected 
roles. The analysis of this paper will be extended to other processes, hke semi- 
inclusive Deep Inelastic Scattering, where many SSA effects have been observed 
[17, 18] and are being measured. 
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Figure 1: Maximised values of vs. xp, at y/s ^ 19.4 GeV and fixed pt = 1-5 
GeV/c, as given by the Collins mechanism alone; the values shown are obtained by 
saturating all bounds on the unknown soft functions and by adding constructively 
all different contributions. The four plots correspond to different choices of the 
distribution and fragmentation functions used to saturate the bounds, as indicated 
in the legends. In each plot the upper, middle and lower sets of curves and data 
refer respectively to n'^, vr^ and tt^. Data are from Ref. [19]. See the text for further 
details. 



21 



